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1 INTRODUCTION 


ABSTRACT 

We report the serendipitous discovery of a late-type giant star that exhibited a smooth, eclipse- 
like drop in flux to a depth of 97 per cent. Minimum flux occurred in April 2012 and the 
total event duration was a few hundred days. Light curves in V, J and K, from the Optical 
Gravitational Lensing Experiment and VISTA Variables in the Via Lactea surveys show a 
remarkably achromatic event. During 17 years of observational coverage of this source only 
one such event was detected. The physical properties of the giant star itself appear somewhat 
unusual, which may ultimately provide a clue towards the nature of the system. By modelling 
the event as an occultation by an object that is elliptical in projection with uniform transparency, 
we place limits on its physical size and velocity. We find that the occultation is unlikely to 
be due to a chance alignment with a foreground object. We consider a number of possible 
candidates for the occulter, which must be optically thick and possess a radius or thickness in 
excess of 0.25 au. None are completely satisfactory matches to all the data. The duration, depth 
and relative achromaticity of the dip mark this out as an exceptionally unusual event, whose 
secret has still not been fully revealed. We find two further candidates in the VVV survey and 
we suggest that these systems, and two previously known examples, may point to a broad class 
of long period eclipsing binaries wherein a giant star is occulted by a circumsecondary disc. 


Key words: stars: peculiar — (stars:) binaries: eclipsing — stars: individual: epsilon Aur — stars: 
individual: TYC 2505-672-1 — stars: individual: ASASSN-21co 


ferometric observations of the 2009-2011 eclipse, including superb 
reconstructed images of the occulting disc in silhouette. 


€ Aurigae is a binary star on a 27 year orbit. The primary is a 
FO type star with supergiant gravity thought to be the result of a 
late thermal flash; the B-type secondary hosts a large and opaque 
disc. The eclipse of the primary by the secondary causes a long 
(^ 2 year), deep (Ay mag ~ 1), flat-bottomed dip in the light curve. 
Kloppenborg et al. (2010, 2015) presented the results of their inter- 
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Denisenko et al. (2013) identified the "optical anti-transient" 
TYC 2505-672-1, and further work by Lipunov et al. (2016) and 
Rodriguez et al. (2016) presented it as a more extreme eclipsing 
binary system than e Aurigae, which had previously held the record 
for longest orbital period. TYC 2505-672-1 is a binary system on a 
69 year orbit, comprising an M-type giant primary with a proposed 
B-type sub-giant companion hosting a large, opaque circumstellar 
disc. The flat-bottomed eclipse lasts — 3.45 years and reaches a 
depth of Ay mag ~ 4.5. 
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A third example of a long-period eclipsing binary analogous to 
€ Aurigae and TYC 2505-672-1 is ASASSN-21co (Way et al. 2021), 
an M-type giant star which began an eclipse event in February 2021. 
A previous event that lasted ~ 80 days was observed as part of the 
All Sky Automated Survey (Pojmanski 1997) in 2009, suggestive 
of a 12 year period. The current event, ongoing at the time of 
preparation of this article, is therefore likely to continue until the 
end of April 2021. 

When it comes to the modelling of occultation events, V1400 
Centauri (also known as Mamajek's Object) springs to mind. It is a 
K5 pre-main sequence star which exhibited a long, deep and com- 
plex dipping event in SuperWASP photometry in 2007. Mamajek 
et al. (2012) found the light curve was consistent with an occultation 
by a large (« 0.4 au) and complex disc and ring system hosted by a 
low mass companion orbiting the primary. 

Another a remarkable occultation event was identified and 
modelled by Rappaport et al. (2019), who also provide a useful list 
of dippers of various types in the literature. It is a single 80 per 
cent depth event of width 1 day in the light curve of the young, 
nearby M dwarf star EPIC 204376071. The dip is notably smooth 
and asymmetric, a trait which the authors ascribe to an occultation 
by an inclined and tilted disc likely hosted by a second body in orbit 
around the M dwarf. 

The VISTA Variables in the Via-Lactea (VVV; Minniti et al. 
2010) survey, and its extension the V VV eXtended (VVVX; Minniti 
2018) survey, sample > 10? sources in the southern Galactic disc 
and bulge over hundreds of epochs across nearly a decade. It is not 
surprising that among them we find many unique variable objects 
exhibiting extreme behaviour that does not fit any known class of 
stellar variability. These special variable sources have been dubbed 
“what-is-this", or *WIT" objects, and follow-up investigations have 
revealed the nature of a few of them. The list of “WIT" objects 
includes: 


VVV-WIT-01, probably an obscured nova or a protostellar colli- 
sion (Lucas et al. 2020); 

VVV-WIT-04, an extragalactic source thought to be an infrared 
violent variable quasar (Saito et al. 2019b); 

VVV-WIT-06, a transient of unknown nature that may be a rare 
Galactic nova, or a stellar merger (Minniti et al. 2017; Banerjee 
et al. 2018); 

and VVV-WIT-07, an intriguing variable source thought to fit 
within the dipper paradigm (Saito et al. 20192). 


Here, we present another unusual transient object — 
VVV-WIT-08 - exhibiting a single smooth, near-symmetrical dip- 
ping event over approximately 200 days. The event is largely achro- 
matic across optical and near-infrared wavelengths. It is at first 
glance reminiscent of a grazing exoplanet transit, but it has a depth 
of ~ 97 per cent in the available wavebands. 

Subsequent shallow searches in the VVV and VVVX data 
revealed two other promising candidate events exhibiting similar 
deep reductions in flux, albeit with poorer light curve sampling, that 
could not be readily explained through common causes of stellar 
variability. Since their data are not as complete as for VV V-WIT-08 
we postpone their analysis and a more comprehensive search for 
other similar events for a future work. 

In Section 2 we summarise the data we have collected 
for VVV-WIT-08. We estimate the physical parameters of 
VVV-WIT-08 in Section 3. Section 4 analyses the event itself and 
describes our attempt to model it, where Section 5 discusses how 
the models compare to the parameters of known objects. We draw 
conclusions in Section 6. 


2 DATA 


VVV-WIT-08 was initially identified during validation of a prelim- 
inary VVV and VVVX light curve database and separately in a 
search for AK, > 4 mag high amplitude variable stars! in that same 
database. 

A further shallow, targeted search in the preliminary database 
yielded a second source at witha AK, « 3 mag event, VV V-WIT-10. 
Subsequent to this, data processing and pipeline improvements led 
to arevised database and a further search for VVV variable stars and 
transients having AK, > 4 mag (Lucas et al., in prep). This further 
search yielded a third source exhibiting a significant unexplained 
drop in flux, VVV-WIT-11. 

The epoch J2016.0 coordinates of WVV-WIT-10 are 
œ = 15h09m33.80s and 6 = -59d38m42.56s, courtesy of Gaia 
eDR3? (Gaia Collaboration et al. 2020). It has Ky = 12.78 mag 
and J — K, = 2.01 mag in the out-of-event photometry. Galactic 
coordinates are /, b = 319.58°, —1.32°, notably a disc field. 

The epoch J2014.0 coordinates of VVV-WIT-11 are 
œ = 17h33m41.26s and 6 = -35d44m16.74s, courtesy of the 
VVV survey. VVV-WIT-11 is present in Gaia eDR3°*, but the 
VVV position is more precise due to the source having a mag- 
nitude of 20.44 in the Gaia G band. It has K, — 12.15 mag and 
J — Ks = 2.38 mag in the out-of-event photometry. Galactic coor- 
dinates are /, b = 352.90°, —1.49*. 

The discovery of VVV-WIT-11 was fortuitous, as the source 
disappears completely for an entire observing season. As a result it 
would not ordinarily have been selected as part of the AK, > 4 mag 
search (which did not include upper limits), were it not for a handful 
of mismatches with an adjacent slightly blended faint star during 
the dimming event. The absence of a detection during the dimming 
event suggests AK, z 5 mag. 

The multi-band time series photometry for VVV-WIT-08, 
VVV-WIT-10, and VVV-WIT-11 are included in the supplemen- 
tary material. Also included is a figure showing the VVV Ks and 
OGLE / band light curve of VVV-WIT- 10. 

A summary of the photometric and other data we have collected 
for VV V-WIT-08 is provided in Table 1 and these data are described 
in more detail below. 


2.1 Light curves 


A detailed description of the VVV and VVVX light curve database 
and data processing is outside the scope of this paper. In summary, 
VVV and VVVX images were processed with a modified version 
of the DoPHOT profile fit image reduction software (Schechter 
et al. 1993; Alonso-García et al. 2012). The resultant photometric 
measurements were then calibrated through a globally optimized 
zero point plus illumination map model, and a further high spatial 
frequency noise reduction routine is used to reduce scatter in the 
light curves. Both photometric calibration stages are described in 
detail in Smith et al. (in prep). VVV and VVVX observations cover 
the 2010-2019 calendar years. 

The Optical Gravitational Lensing Experiment (OGLE) survey 
covers the area of sky containing VV V-WIT-08 as part of OGLE-III 
and -IV (Udalski 2003; Udalski et al. 2015) and provides data from 


! The final measured K, amplitude was a little under 4 mag as the prelim- 
inary AK, > 4 mag measurement included a poor quality observation at 
minimum flux which we've since discounted. 

? VVV-WIT-10 Gaia eDR3 source, id = 5876552066041521024 

3 VVV-WIT-11 Gaia eDR3 source, id = 5974812366541771520 
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Table 1. Data for VVV-WIT-08 obtained from a variety of sources. Pho- 
tometric errors, particularly ground based ones, are likely to be subject to 
additional absolute calibration errors. 


Parameter Value Unit Source 
ID 4044152029396602880 
Q2016.0 18h00m47.41s 
Ó2016.0 -30d36m09.28s 
T 0.0091 + 0.0737 mas 
Haws -3.670 + 0.081 mas yr™! Gaia eDR3 ^ 
HS —6.981 + 0.056 mas yr! 
Gpp 17.3964 + 0.0139 mag 
G 16.0340 + 0.0010 mag 
Grp 14.7494 + 0.0237 mag 
r 16.292 + 0.004 mag 
i 14.974 + 0.003 mag VPHAS+ ^ 
Ha 15.776 + 0.005 mag 
g 18.267 + 0.009 mag 
i 15.354" mag PS1 € 
Z 14.768 + 0.001 mag 
y 14.410 + 0.003 mag 
18.290 + 0.006 mag 
r 16.200 + 0.006 mag 
i 15.145 + 0.006 mag DECaPS 4 
Zz 14.506 + 0.007 mag 
Y 14.234 + 0.006 mag 
Z 14.100 + 0.013 mag 
Y 13.448 + 0.013 mag 
J 12.778 + 0.003 mag VVV e 
H 11.842 x 0.006 mag 
Ks 11.500 + 0.002 mag 
J 12.767 + 0.035 mag 
H 11.736 x 0.031 mag 2MASS Í 
Ks 11.421 + 0.027 mag 
wi 11.118 + 0.031 mag 
2 11.525 + 0.02 
w 323 E 0:028 HUE AIIWISE € 
W3 > 11.668 mag 
W4 > 8.834 mag 
3.6 um 11.292 + 0.048 mag 
4.5 11.255 + 0.064 
uie jene AE GLIMPSE ^ 
5.8um 11.278 + 0.058 mag 
8.0 um 11.193 + 0.045 mag 


a Gaia Collaboration et al. (2016, 2020). 

> Drew et al. (2014). 

© Chambers et al. (2016). 

4 Schlafly et al. (2018). 

* Minniti et al. (2010). 

f Skrutskie et al. (2006). 

8 Wright et al. (2010); Cutri et al. (2013). 

h Spitzer Science (2009). 

! No error given due to this measurement being based on one observation. 
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2001 to 2019. Due to heavy blending and huge brightness variations 
of VVV-WIT-08, the standard pipeline photometry was insufficient 
and a re-reduction and calibration was necessary. V V V-WIT-08 was 
still faintly visible in images taken during the peak of the dimming. 
To produce a high quality light curve at minimum flux, we selected 
good-seeing in-event images and constructed new templates for a 
dedicated image subtraction in the / and V bands. 

The resultant V, J and Ks band light curves for VV V-WIT-08 
are shown in Figure 1. Ks band image cutouts demonstrating the 
high source density of the field are shown in Figure 2. 


2.2 Survey photometry 


In addition to the V and J band optical and Ks band near-infrared 
light curves from OGLE and VVVX respectively, we acquired sur- 
vey photometry from a number of additional sources, as presented 
in Table 1. Note that since the optimal extraction of photometry 
from OGLE and VVVX images required extra care to address the 
nearby sources, so the target is likely to be blended to some degree in 
low angular resolution data such as GLIMPSE, longer wavelength 
WISE bands and possibly 2MASS. 


2.3 Spectroscopy 


Spectroscopic follow-up of VV V-WIT-08 was performed at the 4- 
meter SOAR Telescope in Chile on June 22, 2019, during the Sci- 
ence Verification phase of TripleSpec 4.1 IR spectrograph (Wilson 
et al. 2004; Herter et al. 2008). The instrument covers simultane- 
ously the wavelength range of 0.80 um to 2.47 um with resolving 
power of R~3500 and 1.1” wide slit. The collected spectra were re- 
duced with a modified version of the IDL-based SpexTool (Cushing 
et al. 2004) and flux calibrated against HD 162220. The resulting 
spectrum is presented in Figure 3. 


3 THE NATURE OF VVV-WIT-08 
3.1 Spectroscopic analysis 


A visual comparison of our near-infrared spectrum of V V V-WIT-08 
to those of the IRTF spectral library (Rayner et al. 2009) suggests 
that the star is a cool giant with a type no earlier than K7 and no 
later than M2, due to the general pseudo-continuum shape and the 
lack of FeH at 0.97 um respectively. Additionally, the CN band at 
1.1 um is too weak to be consistent with a supergiant, as is evident 
in Figure 3. 

A separate spectroscopic modelling by interpolated stellar tem- 
plates from Allende Prieto et al. (2018) yielded a best fit effec- 
tive temperature of 3614 + 3 K, [M/H] of —0.29 + 0.01, log g of 
1.05 + 0.01, and a heliocentric radial velocity of -14 + 10 km sl, 
The fitting was performed using a global optimization algorithm, 
implemented in the FERAE code* (Allende Prieto et al. 2006). The 
provided uncertainties are statistical, but systematic uncertainties 
are likely significantly larger due to the relatively high signal-to- 
noise ratio of the near-infrared spectrum. Reasonable estimates for 
the systematic uncertainties based for a typical cool giant are 50 K, 
0.5 and 0.1 for Taf, log g, and [M/H] respectively. 


4 FERAE is available from http://github.com/callendeprieto/ 
ferre 
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Figure 1. OGLE V and / bands and VVV and VVVX Ks band flux relative to the inverse variance weighted mean of the flux outside the event (|t—to | > 365 days) 
for VV V-WIT-08. Inset is a closer view of 1 year of data centered on the minimum 7 band flux measurement. 


Figure 2. VVV K, band image cutouts from three days before minimum 
flux (left) and approximately two years later (right). Cutouts are 10" on 
each side. The position of VV V-WIT-08 is indicated by the circular marker, 
the positions of other sources within 2" are indicated by crosses. VVV, 
DECaPS, and Gaia eDR3 positions are indicated in red, cyan and orange 
respectively. The white space on the left panel is due to the proximity of the 
target to the edge of the detector. 


3.2 Photometric analysis 
3.2.1 Isochrone fitting 


To estimate the physical properties of the VV V-WIT-08, we fit its 
spectral energy distribution (SED) using photometric data from 
various surveys presented in Table 1 and Figure 5. 

We compare the photometric measurements with uncertain- 
ties presented in Table 1 to MIST isochrone predictions (Dot- 
ter 2016; Choi et al. 2016) shifted to various distances and with 
effects of interstellar extinction included. In the modelling pro- 
cedure we also take into account the log g, 7.g and [M/H] con- 
straints from the near-infrared spectrum, and Gaia eDR3 parallax 
measurements. To take into account possible systematic errors in 
isochrones/photometry/reddening we also include in the model an 
additional uncertainty (or excess noise) that is shared across surveys 
and is added in quadrature to each photometric measurement error. 

The full model includes the following parameters: 1) Stellar 
age 2) Metallicity 3) Initial mass 4) Distance 5) Ay extinction and 
6) Excess photometric noise. The isochrone interpolation between 


metallicity, age, and initial mass is performed using MINIMINT? code 
(Koposov 2021). 

To properly estimate the effects of extinction we use the AT- 
LAS9 MONI log g = 1.5 solar metallicity model spectrum (Castelli 
& Kurucz 2003) and the filter transmission from the SVO Filter 
Profile service (Rodrigo et al. 2012) to estimate the effective wave- 
length of each filter for our target. These are then used to estimate 
A4/ Av for the optical and near-infrared photometric bandpasses us- 
ing equations 9 and 10 of Wang & Chen (2019). For the mid-infrared 
data equivalent relations are not provided, so we adopt the Aj/A, 
values for a typical red clump giant provided in their table 3. These 
adopted mid-infrared values should be an adequate approximation 
for our target at these wavelengths due to their relative insensitivity 
to the optical extinction. 

We adopt the following priors on model parameters. The stellar 
age is uniform between 10 Myr and 13 Gyr, metallicity is a Gaussian 
distribution with mean —0.29 and standard deviation 0.1 (based on 
the spectroscopic constraints), extinction Ay is uniform and within 
0 < Ay < 4. The prior on additional photometric scatter across 
different bands is log-uniform between 0.01 and 0.3 mag, and for 
mass we use the prior based on the Chabrier (2005) initial mass 
function with Min = 0.5Mo and Mmax = 10Mo. 

To sample the posterior of the model described above we use 
parallel tempering ensemble sampling, implemented by the ptemcee 
package (Vousden et al. 2016; Foreman-Mackey et al. 201 3)6, with 
1000 walkers and ladder of 5 temperatures that were run for 50000 
steps after a burn-in of 50000 steps which were discarded. 

The preliminary isochrone modelling outcome, shown in Fig- 
ure 4, suggests that the giant can be described by 1.0 — 2.4 solar 
mass star at a range of distances from ~ 17 to 42 kpc at 2 e. How- 
ever, an important additional consideration is the range of allowable 
tangential velocities. It is extremely unlikely that VV V-WIT-08 has 
a space velocity in excess of the Galactic escape velocity. If we 
require that Vian(R) < Vesc(R) then we can constrain the giant pa- 
rameters further. Vesc(R) is obtained using the AGAMA software 
package (Vasiliev 2019) with the McMillan (2017) potential and the 
distance (R) samples from the modelling. Vesc is between ~ 500 
and ~ 600 km s^! at these distances and in this direction. We ignore 


5 https://github.com/segasai/minimint 
6 nttps://github.com/willvousden/ptemcee 
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Figure 3. A TripleSpec near-infrared spectrum of VVV-WIT-08 obtained on the 22 June 2019 (dark blue line). The spectrum is smoothed with a narrow 
(^ 6.5 A) Gaussian filter for clarity. In addition we include four template spectra of stars of a similar Ty but different luminosity classes, obtained from the 
IRTF Spectral Library (Rayner et al. 2009). The templates were normalised to the mean flux of VVV-WIT-08 in the Ks band, reddened using Ay = 2.17 mag 


(see Section 3.2 for details), and then offset for clarity. 


Table 2. Priors and median and 95th percentile credible intervals on the 
posterior pdf of the physical properties of VVV-WIT-08 on requiring that 
Vian (R) < Vesc (R). 


Parameter Prior Posterior Units 
[M/H] G(-0.29, 0.1) —0.18*0-16 dex 
-2 4.37 
Age U(1074, 13.0) 8.34* 06 Gyr 
Mass I MF(0.5, 10.0)@ (40 Mo 
0.28 
Av Uu (0, 4) 1.677 056 mag 
Distance LU(1, 100) 19.9*L7 kpc 
c I4 (1072,03) 0.077005 — mag 
Test 3659+75 K 
0.18 
log g 0.8175 33 dex 
w ag. nas 
i , 0.2 3 
Luminosity 0.7755 10° Lo 
; 0.03 
Radius (rg) 0. 3275. e au 
Radius (rg) 15. P i nas 
VI —493+140 — kms! 
Vp -15* — kms! 


@ Chabrier (2005) IMF prior 


the relatively small component of radial velocity for this estimate. 
Figure 4 shows the impact of this requirement on the posterior prob- 
ability distribution, and Table 2 summarises the inferred parameters 
and derived quantities. Essentially the part of the posterior with he- 
liocentric distances larger than ~ 20 kpc is removed. These results 
suggest that the giant is likely a 0.9 — 1.4 solar mass star at a he- 
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Figure 4. Posterior pdf on giant physical parameters obtained by modelling 
the photometry with MIST isochrones. The blue distribution is allowing an 
unrestricted tangential velocity, the green requires that the tangential velocity 
is lower than the Galactic escape velocity. 


liocentric distance of between 16.1 and 21.5 kpc. The photometry 
of the model matches well the observations across all the bands 
requiring an additional systematic error of only ~ 0.07 mag. We 
adopt the Vian(R) « Vesc(R) requirement for all further analyses. 
Inthe upper panel of Figure 5, we show the best fit sample from 
the isochrone spectral energy distribution (SED) model versus the 
survey photometry. Also shown are the 68 and 95 per cent spread 


6 L.C. Smith et al. 


O 2MASS o DECam o Gaiao IRACO PSo VISTA oO WISE 


o 
H 
1 
poe 
-e 
et 
re 
æ 
Heie 
me 
m 
—— 
Ea 


g 0.04 | | 
al TTT dH | i 
i i | 
-0.2 4 
0.5 0.6 0.8 1 2 3 4 5 6 7 8 
Aerr | um 


Figure 5. Upper: Isochrone SEDs of a low y? sample from the posterior 
of the physical parameters of VV V-WIT-08. The black line shows a low x? 
model. Lower: The offset from the median and spread in synthetic photom- 
etry, 95 and 68 per cent confidence intervals and the median are shown. The 


modelled excess error component o = 0.07:0.05 mag is not included. 


in the offsets relative to the photometric data points measured from 
a selection of random samples. 

Crucially, from the above analysis we are able to obtain an 
estimate of the luminosity and hence radius of the giant. Evaluat- 
ing the isochrone interpolant at each posterior sample, we obtain 
the median and 95 per cent confidence interval on luminosity of 
rs x10? Lo. From the luminosity and 7. samples, we obtain 
a median and 95 per cent confidence interval on the stellar radius 


+0.03 
of 0.3275 o; au. 


3.2.2 Radius at Galactic bulge distance 


Taking the Gaia eDR3 proper motion of the source and the 
isochrone fitting distance to the giant (with the Vtan(R) < Vesc(R) 
requirement), we obtain tangential velocities of —493+140 and 


-15+ km s^! in the Galactic longitude and latitude directions 
respectively, having corrected for solar reflex motion. The lower 
bound on V; brings it closer to that of the Galactic rotation curve 
but it's still larger than might be expected. The distance distribution 
of stars from a Galaxia (Sharma et al. 2011) model at the sky 
position of VV V-WIT-08 that also share its proper motion within 
its uncertainties is shown in Figure 6. Note that a distance below 
14.2 kpc is ruled out to 3 c from the isochrone analysis in Section 
3.2.1 under the assumption that VVV-WIT-08 is on an ordinary 
stellar evolutionary pathway. This is marginally consistent with the 
inference based on Galaxia. 

The space velocity of VV V-WIT-08 implied by the analysis in 
Section 3.2.1 appears halo-like, though this is in tension with its 
relatively metal rich spectrum ([M/H]- —0.29). We note that Vp 
and Vos are more moderate. These velocity components and the 
metallicity appear more consistent with a disc or bulge membership. 

Since the proper motion of the star from VVV and VVVX 
astrometry based on Gaia DR2 stars as astrometric reference points 
matched well the one measured by Gaia within 20 and therefore is 
likely accurate, we consider that the physical parameters of the giant 
star obtained through fitting of the isochrones to the photometry 
are relatively unlikely prima facie. We note that in addition to the 
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Figure 6. The pdf of the various independent and dependent parameters 
of the bulge distance physical parameter analysis. The distributions of Teg 
and Ay are from the analysis of the spectrum of VV V-WIT-08 and nearby 
red clump stars respectively. The distance distribution is that of stars in a 
Galaxia model with similar position and proper motion to VVV-WIT-08. 
The Ks band measurement assumes 0.02 mag absolute calibration error. 
The distributions of luminosity and radius are extracted from 107 random 
samples from the other parameters and model giant spectra. 


metallicity, the proper motion is consistent with those of Galactic 
bulge stars. 

To assess what distance would be expected from this proper 
motion alone we use the Galaxia simulation at the coordinates of 
VVV-WIT-08 described above, having selected stars that share its 
proper motion on taking into account the Gaia eDR3 uncertainties. 
Other than a tiny component of nearby disc stars, which we remove, 
there are two modes in the distance distribution corresponding to 
stars in the bulge and the far side of the disc. The heliocentric 
distance range of 6.0 to 13.3 kpc contains 95 per cent of the Galaxia 
sample. 

The difference between the kinematic and isochrone dis- 
tances may be interpreted as as lower-than-expected luminosity of 
VVV-WIT-08 than predicted by isochrones and can perhaps be ex- 
plained if we postulate that it did not follow a typical evolutionary 
pathway, as a fit to the MIST isochrones would require. There are 
examples in the literature of under-luminous giant stars, namely 
sub-subgiant stars or SSGs (Geller et al. 2017a; Leiner et al. 2017; 
Geller et al. 2017b). These are redder than ordinary main sequence 
stars and fainter than giant and subgiant stars. SSGs are thought to 
form through binary evolution channels, although other formation 
mechanisms have been proposed (Geller et al. 2017b). The majority 
of SSGs are variable to a significant degree, either in their pho- 
tometry or radial velocity. While VV V-WIT-08 is not significantly 
variable in the time series photometry we have (see Section 3.3) — 
barring the dipping event itself — radial velocity follow up would 
be worthwhile. Approximately one third of SSGs are Ho emitters 
(Geller et al. 20172), but this does not appear to be the case for 
VVV-WIT-08 in the VPHAS+ photometry. 
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The modelling of exotic stellar evolutionary channels is dif- 
ficult. An attempt to find suitable parameters using BPASS binary 
evolution models (Stanway & Eldridge 2018) could not be recon- 
ciled with the ordinary atmospheric parameters obtained from the 
spectrum of VV V-WIT-08 and hence was fruitless. With this caveat 
in mind, we attempted to estimate the radius of VVV-WIT-08 by 
considering simply the luminosity required to obtain the measured 
Ks band magnitude given the Tz, log g, and [M/H] parameters ob- 
tained from the NIR spectrum and the distance pdf obtained from 
the Galaxia model as described above. We include an additional 
0.02 mag uncertainty to take into account the typical absolute cali- 
bration error. The Ks flux measurement should be the least impacted 
by uncertain extinction without resorting to using the mid-infrared 
photometry which is likely to suffer from blending. Through a mea- 
surement of E(J — Ks) = 0.36 +0.04 mag for nearby red clump stars 
using VISTA PSF photometry and the Wang & Chen (2019) law we 
obtain an estimate of Ay — 2.17 x 0.26 mag for this field, which 
we adopt. We used appropriate ATLAS model spectra (Castelli & 
Kurucz 2003) and the synphot python package (STScI Develop- 
ment Team 2018) to perform synthetic observations and estimate 
stellar luminosities and hence radii. Sampling from the probability 
distributions of these parameters 107 times yields a radius distribu- 
tion of 0. 14+0-07 au Gi Ro, median and 95 per cent confidence 
interval). Figure 6 shows the probability density functions of all 
parameters. 

The measurements made assuming a Galactic bulge location 
will be retained in addition to the isochronal measurements for the 
purposes of further analysis. We consider the former to be sensible 
lower limits on the distance and (more importantly) radius of the 
giant, where the latter are sensible upper limits. It’s highly likely that 
the distance and radius of VVV-WIT-08 lies somewhere in between 
these two scenarios. 

We note that if we adopt a Galactic bulge location, an alterna- 
tive explanation for the low luminosity of the giant would be that 
there is diffuse matter in the star system or along the line of sight 
that is causing grey extinction. We have no reason to think this is 
correct but given the potential presence of optically thick material 
causing a large drop in flux it cannot be ruled out. In this scenario 
it would be appropriate to adopt the isochronal radius distribution, 
which is the most important parameter for the remainder of this 
investigation. 


3.3 Quiescent variability 


Close inspection of the light curve presented in Figure 1 reveals a 
quasi-periodic signal in the / band data outside the event. We also 
find it is necessary to inflate the OGLE 7 and V band photometric 
uncertainties through the addition in quadrature of an excess error 
Cexcess,V,1 = 0.019,0.012 mag in order to explain the scatter in 
their extra-event light curves. Uncertainties in the Ks band data 
for this source appear to be overestimated by approximately 40 per 
cent at baseline flux. We suspect DoPHOT is overestimating the 
uncertainties near saturation that occurs close to Ks = 11.5 mag in 
VVV and VVVX observations with good seeing. 

It remains possible that this apparent variability could be ex- 
plained by residual systematic errors in the data, but we consider 
this less likely given the care taken when reducing these data. We 
also note that the V and / band scatter is mildly correlated. 

To attempt to measure a period in the variability while taking 
into account possible time-varying stellar surface conditions, we 
produce mean Lomb-Scargle periodograms across moving windows 
in the extra-event J band data using a sequence of different window 
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widths. We find a weak period of approximately 195 days (and its 
alias at approximately 390 days) and a weaker one at approximately 
124 days. If genuine, these periods would be roughly consistent with 
rotation speeds of red giants measured with Kepler, in particular a 
124 day period would fit nicely within the distribution shown in 
figure 5 of Ceillier et al. (2017). The magnitude of the excess scatter 
appears to be consistent with that of ordinary variability of some red 
giant stars as measured by Henry et al. (2000), and Oláh et al. (2018) 
found that approximately half of giants exhibit flux variability above 
1 per cent. 

A ‘back of the envelope’ estimate of the flux variability ex- 
pected from a single star spot gives around 1.8 per cent in the V 
band and 1.3 per cent in the 7 band, assuming black body spectra 
and requiring AT = 200 K with a spot covering 5 per cent of the 
giants surface. For reference, the variability in the Ks band using 
this model is 0.6 per cent. 

It is worth bearing in mind that the above tentative period of the 
quiescent variability is relatively close to the duration of the major 
dimming event. This should be taken under consideration during 
analysis of the event itself. In particular we suspect that the subtle 
difference in the flux at baseline immediately before and after the 
event is due to this quiescent variability. 


4 EVENT ANALYSIS 
4.1 Intrinsic variability 


We are not aware of any mechanism of intrinsic stellar variability 
which could cause the observed drop in flux of this magnitude across 
the wavelength range 0.5 — 2.2m. The shape of the dipping event 
and relative constancy of its colour are additional strikes against 
the intrinsic variability hypothesis. We briefly discuss two specific 
possibilities, before discarding them. 


4.1.1 R Coronae Borealis-like variability 


R Coronae Borealis (RCB hereafter) are rare, hydrogen-deficient 
carbon supergiant stars. Their light curves are characterised by deep 
(up to 9 mags in V band) reductions in their luminosity due to the 
condensation of carbon to soot in their atmospheres. RCB stars have 
warm shells with 300 < Tag (K) < 1000 and hence tend to exhibit a 
mid-infrared excess (Tisserand 2012), whereas this does not appear 
to be the true of VV V-WIT-08. Major dimming events of RCB stars 
are characterised by a sharp ingress and a relatively shallow egress. 
In principle, the VV V-WIT-08 event light curve is consistent with 
this due to the gap in the ingress data but its relative achromaticity 
is not - RCB dimming events are significantly deeper at shorter 
wavelengths (Feast et al. 1997). 


4.1.2 Variable young stellar object 


High amplitude photometric variability is common in YSOs, often 
caused by episodic accretion. Dips are also known to occur. Low 
amplitude dips typically have timescales from hours to days (Find- 
eisen et al. 2013; Cody et al. 2014; Scaringi et al. 2016; Ansdell 
et al. 2016) but deeper dips on timescales of years are also seen 
(e.g. Bouvier et al. 2013; Contreras Peña et al. 2017), attributed to 
a substantial perturbation of the accretion disc structure, perhaps 
by a young planet. However, VV V-WIT-08 does not appear to re- 
side in or near a region of the sky with obvious signs of ongoing 
star formation. Its spectrum shows no signs of youth, and its lack 
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of significant mid-infrared excess indicates similarly. Outside the 
event, VVV-WIT-08 has low variability and the dimming event it- 
self is largely achromatic, both of which are uncharacteristic of YSO 
variability. 


4.2 Occultation 


As we are unable to explain the behaviour of VV V-WIT-08 through 
intrinsic variability, we turn our attention to occultation. 


4.2.1 Preliminary considerations 


If an occulter is orbiting the giant, then the long OGLE / band time 
baseline rules out an orbital period shorter than around 11 years. 

The time baseline between the first point of reduced flux and 
the approximate last point of reduced flux, and hence the minimum 
event width, is around 170 days. Additionally, in this sight line 
and due to the distance to the giant, some portion of the Galactic 
bulge is in the foreground. The long timescale of the event, and the 
possibility that the occulter is at an intermediate distance means that 
relative parallax between the giant star and the occulter may need 
to be considered. 

As might be expected given its position and as is demonstrated 
in Figure 2, the target is subject to significant crowding. The nearest 
on-sky objects are 1.06” and 1.44" away in the Gaia DR2 data, and 
there is an additional source at 1.74" separation in the VVV and 
DECaPS data which was not detected by Gaia. None of these neigh- 
bours exhibit clear dimming events similar to that of VVV-WIT-08 
in our data. This means any occulter must either be more distant 
than the nearby objects, or have angular size in the direction of the 
nearby stars that is smaller than their separations and a motion vec- 
tor which does not pass in front of them during our observational 
coverage. 

The lack of obvious amplification in the light curve due to 
microlensing implies that the angular size of an occulter must be 
much larger than its Einstein radius. 

There is a several hundred day gap in our observational cov- 
erage during most of the 2009 observing season. In principle, a 
similar event could have occurred in this year and been missed, but 
a periodic signal of « 3 years is ruled out due to the lack of similar 
events in 2003, 2006, and 2015. 

Due to surface brightness considerations, the depth of the event 
rules out a common eclipsing binary system with the second com- 
ponent being another star. To produce a 97 per cent reduction in 
the / band flux the secondary would need to have an equally large 
radius and an effective temperature below 1400 K. The only remain- 
ing binary system parameter space is a large (r z rg) and dark or 
self obscured companion. 


4.2.2 Chance alignment probability 


If the occultation is due to the presence of an object unbound from 
the background star at some intermediate distance then we are able 
to estimate the implied space density required in order to have a 
reasonable chance of our finding VVV-WIT-08. On adopting the 
approximation that occulters and background stars are circular in 
projection and ignoring the effects of parallax motion, if we have £o 
surface density of occulters and Xs surface density of background 
stars, moving with relative velocity V with respect to occulters, then 
the rate of encounters with minimum separation between r and r+dr 


per unit time dt occurring in the area A will be 
dR = 2VXo%Xs5 Adrdt. 


Occultation duration and depth are functions of V and the radii of the 
occulter and background star (ro, rs). We consider depth purely as 
a fraction of the background stellar disc that is obscured, ignoring 
limb darkening for simplicity, and ignoring the potential for flux 
contamination or a semi-opaque occulter. 

An estimation of rates requires knowledge of both the stellar 
and occulter density and the relative velocity distribution, and for 
these we employ a Galactic population model. We ran a Galaxia 
simulation at the location of VV V-WIT-08, selecting objects with 
no regard to apparent or absolute magnitude as potential occulters, 
and those with Ks < 17 (which is the approximate limit of our VVV 
psf data) as potential background sources. Using these selections, 
and the further requirement that the background star is more distant 
than the occulter, we evaluate the rates of occultations of different 
durations (p) and depths (Af; depth of one corresponds to full 
occultation). Each sample uses the Galaxia distances and velocities 
for both components and for the radius of the background star, but 
the occulter radius we fix at 0.5 au. 

The rates obtained for p > 10 day occultations are 0.80 and 
0.73 events per square degree per year for Af > 0.5 and Af > 0.9, 
respectively. For p > 100 day events these rates drop to 1.3 x 107? 
and 1.2x 107? deg ?yr-! for the same A f limits. We note that there 
is significant volume in the distribution at Af = 1.0, which is to be 
expected where ro > rs. If we could be certain that the remaining 
3 per cent flux seen in the VVV-WIT-08 event were due solely to 
an incomplete obscuration then we might additionally incorporate 
the requirement that Af < 1.0. In this case, for p > 100 day events 
with 0.9 < Af < 1.0 the rate is 3.4 x 107^ deg ?yr- |. 

We consider the following additional modifiers to the above 
rates. The Galaxia simulation predicts a Ks < 17 mag source den- 
sity of 106 stars per square degree and the rates scale proportionally 
relative to this. The rates reported are under the assumption that the 
occulter space density is equal to the stellar space density. Changing 
the adopted occulter space density will alter the rates proportionally. 
Provided that rg > rs, the rate of all occultations can be directly 
scaled up or down to take into account larger or smaller occulters, 
although the depth and duration distributions will change. 0.5 au is 
appropriate for VV V-WIT-08 (see Table 4). 

Ignoring the b « —5 deg region, where the source density is 
much lower, the VVV survey bulge component is ~ 200 deg ? and 
covers 10 years of observation. These 10 years of observation are 
punctuated by ~ 6 month seasonal gaps in observation, reducing 
detection efficiency by ~ 50 per cent. In our VVV psf data we 
find a source density that is approximately 50 per cent larger than 
Galaxia's prediction of 106 deg ?. Using the 1.2x 107? deg ?yr! 
rate from above for p > 100 day, Af > 0.9 events, we should 
expect to find 18 VVV-WIT-08-like events in our data, with the 
important caveat that this assumes the space density of occulters 
is equal to stars. An additional caveat is that the rates above are 
appropriate for background stars that are brighter than K, — 17 mag, 
the approximate detection limit of our data, but stars must in fact be 
brighter than this at baseline such that we can identify a reduction in 
brightness. The exact baseline magnitude at which we are sensitive 
to dips is difficult to determine, but is likely to be a magnitude or 
two brighter than the K, — 17 mag detection limit of the survey. As 
a result the expected event rate will be a factor of a few lower, and 
hence in order to expect to find VV V-WIT-08 the required occulter 
space density must to first order be similar to that of the stellar 
content. 
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While in principle it is possible that a population of dark oc- 
culting bodies exists that is similar in number to the stellar content 
— if so it would take a concerted search to find them (that to our 
knowledge has not been undertaken in earnest), as we have just 
demonstrated — we consider it unlikely. Furthermore, V V V-WIT-08 
is something of an anomaly itself (in either velocity or radius, see 
Sections 3.2.1 and 3.2.2), and a rare event occurring to a rare object 
compounds the unlikeliness. As such, we deem it improbable that 
a chance alignment of an unrelated object is the cause of the dim- 
ming event. We adopt an occulter bound to the giant as our primary 
hypothesis. 


4.2.3 Modelling 


A plausible explanation of a smooth, achromatic dip in flux is an 
occultation by an optically thick body. We first modelled the light 
curve as a solid spherical occulter, passing on a straight line relative 
to the background star. This modelling was implemented using the 
publicly available BATMAN stellar transit light curve modelling soft- 
ware (Kreidberg 2015) and the pyNesty dynamic nested sampling 
package (Higson et al. 2019; Speagle 2020), and yielded a mediocre 
fit to the light curve with ro/rg ~ 1.3. The best fit we were able 
to achieve is shown in Figure 8 (red line), and demonstrates that 
a spherical occulter on a straight line trajectory cannot adequately 
match the data at hand. 

This led us to enhance the model by allowing the occulter to 
be elliptical with finite transparency and to have a parallax with 
respect to the giant star. This is partly motivated by Rappaport et al. 
(2019), who showed that an occultation by an object that is elliptical 
in projection is able to explain a smooth, asymmetric light curve in 
the case of EPIC 204376071. 

We implemented the light curve generation ourselves using 
custom written C++ software that we release with the paper’. 
Specifically, this software integrates over the obscured part of the 
disc of the star, taking into account the limb darkening. The inte- 
gration is done using a pre-calculated grid of 1000x1000 points on 
the disc of the star. The C++ software is optimized for maximum 
performance and is parallelized using OPENMP. 

Analysis of the space density of dark occulters required in order 
to produce a VV V-WIT-08-like event (Section 4.2.2) indicates that 
there is no need to include parallactic motion in the model. Despite 
this, we were curious to see whether there was volume in the pos- 
terior for high parallax models. Considerable effort was expended 
in trying to extract the posterior pdf of the model parameters with 
parallax included. We attempted the modelling with a large selec- 
tion of currently available nested and ensemble sampling packages. 
Despite intensive effort, careful reparametrisation to reduce degen- 
eracies, and in some cases modification of the sampling packages 
themselves, none were capable of satisfactorily capturing the com- 
plex shape of the posterior and usually missed one or more known 
modes. This was principally due to the narrow likelihood surface, 
multiple clearly distinct modes, and vast regions of the posterior 
with a flat likelihood. 

While we know models containing a large relative occulter 
parallax exist that are capable of explaining our data, the analysis in 
Section 4.2.2 indicates these are unlikely a priori. Due to the need 
to further reduce degeneracies in order to produce a satisfactory 
posterior we ultimately chose to remove this parameter from our 
model. In doing so we reduce the motion vector to a straight line 


7 https://github.com/segasai/vvv dipper. code 
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relative to the giant star, which can be parametrised as a scalar 
value, reducing the velocity posterior distribution to a single mode 
and removing the velocity direction parameter. 

The final model that we implement has the following parame- 
ters and priors: 


(i) The semimajor axis rg of the elliptical occulter. It is measured 
with respect to the angular radius of the star rg. The range over 
which we sample is —0.5 < logig rmaj < 2 with the prior Es We 
discuss the need for this prior below. 

(ii) The semiminor axis Fmin of the elliptical occulter. Again mea- 
sured with respect to the angular radius of the star rg. The prior 
range is —0.5 < logj9/min < 2 with an informative Beta prior on 
the axis ratio min/rmaj = 4 ~ B(1,0.5). 

(iii) The orientation angle of the occulter major axis 0, which is 
measured with respect to the velocity vector. We use the prior 0 — 
(4 (0, x). 

(iv) The impact parameter (radius) rg of the occulter relative to the 
giant at time fo, i.e. the minimum separation between the occulter 
and giant measured in units of the radius of the star rg. The prior is 
logio ro ~ U(-3, 2). 

(v) The tangential proper motion vector of the occulter relative to 
the giant, u. The proper motion is measured in units of radius of the 
giant per day u ~ U (0.0, 0.2). 

(vi) The quadratic limb darkening parameters of the giant a and 
b, defined so that the specific intensity on the surface of the star 
I(u) = 1 — a(1— u) - b(1— u)?, where u = cos O and @ is the 
angle between the line of sight and the normal to the star surface. 
a ~ 4(0.5,0.7), b ~ U(0.05, 0.25). Bounds on the quadratic limb 
darkening parameters a and b cover the range of values expected 
for stars with similar Tag, log g, and [M/H] to VV V-WIT-08 in the 
I band as computed by Claret & Bloemen (2011). 

(vii) The time fg at which the position of the occulter is xo, yo. It 
is the time of the minimum on-sky separation between the occulter 
and giant. The prior is tg ~ 14 (54500, 57500) 

(viii) The transparency of the occulter T, where values of 0 or 1 cor- 
respond to a completely opaque or transparent occulter respectively 
and the prior is logjg T ~ U(-5, 0). 


The model parameters, priors and posteriors from the modelling are 
also summarized in Table 3. 

Two aspects of the prior require a special remarks. The rea- 
son for the Que prior was to correct for the fact that large oc- 
culters are intrinsically more likely to produce occultation events 
due to their large cross-section. With the E prior our poste- 
rior is characterizing the information we gained on the occultor 
conditional on the event taking place. The prior on the axis ratio 
Tmin/"maj = 4 ~ B(1,0.5) was introduced to deal with a degen- 
eracy that our data is mostly constraining the semiminor axis of 
the occulter with essentially no constraints on raj, thus allowing 
arbitrarily elongated objects. The prior on the axis ratio penalizes 
such solutions. This prior has ~ 31 per cent of prior volume in the 
q > 0.9 region and 39 per cent prior volume for 0.5 < q < 0.9. 
Since the prior is informative and is affecting the posterior, caution 
must be exercised when interpreting the inference of raj and rin. 

We parametrise the occulter as having transparency between 
negligible and total. Note however that transparency is wholly de- 
generate with flux contamination, the contaminant being the occul- 
ter or an unrelated background or foreground star. In the absence 
of other information, we are unable to differentiate between flux 
contamination and occulter transparency. 

The model specified above could in principle be applied to a 
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Table 3. Model parameters of the eclipse with priors and posterior measure- 
ments. Posterior column contains the median and 95th percentile credible 
interval in the parameters. Positions and proper motions are relative to those 
of the giant star, and are in units of Tg. 


Parameter Prior Posterior Units 
logi Fmaj U(-0.5, 2) 4 0.24022 Tg 
logio Tmin U(-0.5, 2) 0.004*0-010 Tg 

0 U (0, 7) D radians 
logio ro 14 (-3,2) 1.5108 Tg 

u U (0.0, 0.2) 0.01920 00 1.— re day! 
a 4 (0.5, 0.7) 0.65720 

b «4 (0.05, 0.25) 0.187097 

to U (54500, 57500) 56021.2*50 MJD 
logio T U(-5, 0) -1.52 907 


^ Also includes additional terms rt and rgin/rmaj ~ (1, 0.5), see text. 


aj 


data set with photometric measurements in multiple bands while 
requiring different limb darkening parameters (a, b) for each but 
since we have an order of magnitude more data points in the J band 
than in the V and Ks bands, we chose to model the / band data only. 
The model makes a prediction for the fraction of flux of the star as a 
function of time f(t). We then use those predictions to compare with 
the observed fraction of the star's flux with respect to the baseline 
value of J = 14.713 mag. To deal with the fact that the light curve of 
the star outside the event shows a slight flux variability (caused by 
either uncorrected systematics or low-level intrinsic variability, see 
Section 3.3), we inflate the magnitude uncertainties by 0.012 mag: 
Pini = Jod + 0.0122. 

The posterior on the 9 parameters of the occultation model was 
sampled using nested sampling (Skilling 2004) via the pyMultinest 
package (Feroz et al. 2009; Buchner et al. 2014) with 5000 live 
points. The resulting posterior is shown on Figure 7. We see that 
the inference is robust, with two distinct posterior modes corre- 
sponding to a completely opaque occulter, and an occulter with 
a 3 per cent transparency. The opaque occulter requires that it is 
marginally smaller along one axis than the giant star, and implies 
that the remaining 3 per cent flux in the light curve is due to a 
small un-obscured portion of the giant. The semi-transparent occul- 
ter requires that it is marginally larger along one axis than the giant, 
and implies that the remaining 3 per cent flux in the light curve 
is due wholly to the transparency of the occulter. The two distinct 
transparency regimes result in two non-overlapping occulter orien- 
tation (0) modes. Generally speaking, the semi-transparent occulter 
allows more freedom across the remaining parameters. A selection 
of model light curves sampled from the posterior is shown in Figure 
8. The maximum likelihood models from each of the distinct modes 
are illustrated in Figure 9. 


4.3 Event chromaticity 


With a model of the event in hand we are now better able to quanti- 
tatively assess the chromatic nature of the event. 
Obtaining models derived from the J band data but which are 


appropriate for the V and Ks bands simply requires substituting the 
limb darkening parameters, while keeping the geometric parameters 
constant. If our assumptions that the event is achromatic and the 7 
band model are correct, then we should find that the V and Ks 
band data fit these modified models. By interpolating Claret & 
Bloemen (2011) limb darkening parameters we obtain V and Ks 
band equivalents to those of the 7 band models, then we run these 
and the geometric parameters through the occultation modelling 
software described in Section 4.2.3. The resultant distributions of V 
and Ks band models are shown in Figure 10 against the observations. 

The V band data generally agrees with the V band model pre- 
diction within the photometric errors and model spread, indicating 
no significantly different from the model V — J colour. The same 
cannot be said for the Ks band data, which is a few sigma below 
the prediction of the model between the start of the event and its 
peak. When the Ks band observations resume towards the end of 
the event they are back where the model predicts, indicating a grey 
I — Ks colour at this time. The 20 to 30 per cent lower Ks band 
flux around the event peak evident in the lower panel of Figure 10 
and is computed relative to the predicted model flux. If one instead 
considers the residual as a percentage of the baseline flux then it is 
^ ] per cent. 

We remark that the one further event candidate for which we 
have in-event data, VV V-WIT-10, also appears to exhibit a blue 
I — Ks colour during the event. It is not detected in the OGLE V 
band observations. 

One possible explanation for the blue 7 — Ks colour might 
be that flux contamination at baseline by a fainter, hotter object 
is masking the true obscuration fraction. A hot companion would 
contribute more flux in optical passbands, which might then not 
be obscured during the event, leading to an apparently shallower 
eclipse at shorter wavelengths. To test this hypothesis we reran the 
isochrone modelling described in Section 3.2.1, this time included 
contamination from a fainter, hotter blackbody companion. The 
contaminating source is parametrised as a blackbody with a log- 
uniform 5000 < T,g/K < 50000 prior, and a contaminating flux 
fraction in the Ks band (relative to the giant baseline flux) with 
the uniform prior between 0 and 0.1. The experiment revealed no 
clear evidence for such a contaminating component, although it is 
not ruled out with a high degree of confidence. It is unfortunate 
we do not have sufficient data to conduct this test with the in- 
event photometry. With the baseline photometry we are trying to 
identify a signal which is approximately the same size as the typical 
photometric uncertainties. 

Weak evidence in favour of the flux contamination scenario 
is the presence of the 3 per cent transparency mode in the event 
modelling posterior, which is directly interchangeable with flux 
contamination. That is, the measured 3 per cent transparency might 
instead be 3 per cent flux contamination in the 7 band. Further evi- 
dence for the flux contamination scenario is that the colour change is 
expected to be proportional to the fraction of obscured giant flux. It 
therefore is able to explain the return to a zero (or near-zero) I — Ks 
colour at egress. Note that the contaminant need not be the occulter 
itself, it could be another member of the system or an unrelated 
object along the line of sight. 

A second, discounted, hypothesis was that an opacity hole in 
the center of the occulting disc would allow light in from the central 
region of the disc of the background giant. In our occultation model 
the center of the occulting disc has a relatively low impact parameter 
in some cases. As flux is concentrated more towards the limb of the 
background giant in the near-infrared this would manifest as a deeper 
eclipse at these wavelengths. The counterargument to this theory is 
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Figure 7. Posterior pdf on elliptical occultation model parameters. Length parameters (occulter size, position and motion) are in units of giant radii (rg) and 
motions are relative to the giant. Times are in units of modified Julian days and @ is in units of radians. The contours correspond to lo, 20 and 30 limits. 
The median and 95 per cent confidence interval on these parameters are given in Table 3 


that the opposite effect is expected to be seen at shorter wavelengths, 
ie. where the Ks eclipse would be deeper than the 7 band, the V 
ought to be shallower, yet this is not the case. Furthermore, the 
discontinuity in J — Ks colour between the early and late eclipse 
cannot be explained through this method. 


Another possible explanation for the blue / — Ks colour may be 
that some scattered light from the edges of the occulter was detected 
at mid-eclipse. The single scattering albedo of sub-micrometre dust 
grains is higher in the optical waveband than in the infrared, see 
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e.g. Lucas & Roche (1998). The phase function of such grains also 
causes an increase in forward scattering in the optical, relative to the 
infrared. Therefore, if the achromatic nature (to first order) of the 
event is due mainly to high optical depth rather than a large grain 
size distribution, then forward scattering of light from the giant star 
by the edges of an occulting disc can plausibly explain the slightly 
greater dimming in the infrared. This explanation is invalid if the 
disc is transparent to some degree, and while the modelling tends 
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Figure 8. Upper: A random selection of 200 light curve models from the posterior pdf with the OGLE J band measurements from which they were obtained. 
Lower: Fractional residuals to the maximum likelihood estimate (MLE) model. The maximum event depth in the 7 band as measured from these models is 
96:03:07 percent at 1 o. In both panels the light curve models are rendered in blue where they correspond to the high transparency mode, and in green where 
they correspond to the low transparency mode. The error bars do not reflect the additional 0.012 mag uncertainty added in quadrature to take into account 
quiescent variability. The lower panel includes in red the difference between the elliptical occulter maximum likelihood model and the best fit we were able to 


achieve with an occulter that is opaque and circular in projection. 
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Figure 9. An illustration of the elliptical occulter eclipse models using maximum likelihood estimate (MLE) parameters for an opaque occulter (logi; T < —2, 
green) and a semi-transparent occulter (log;g T > —2, blue). 
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Figure 10. Fractional residuals to the median model flux for the V, J and Ks bands. The V and Ks models were obtained by substituting appropriate V and Ks 
band limb darkening parameters into the models obtained by fitting using the / band data. The dark and light blue bands indicate 68th and 95th (respectively) 


percentile spread between the model samples. 


to prefer a 3 per cent transparency, this is not a given (as discussed 
in Section 5.1). 


5 DISCUSSION 
5.1 Event modelling 


Figure 11 shows the posterior pdf on occulter parameters implied 
by the event modelling for the two giant star models considered: the 
isochrone-based parameters with the Vian(R) < Vesc(R) require- 
ment (Section 3.2.1, ‘isochrone model’ hereafter), and the parame- 
ters implied by requiring the giant is located in the bulge (Section 
3.2.2, ‘bulge model’ hereafter). Further details are provided in Table 
4. 

The relative velocity obtained from the modelling of our oc- 
cultation is well constrained by the data and strongly covariant with 
the radius of the giant star. If we assume a circular orbit of the 
occulter around the giant then the only remaining free parameter in 
the determination of orbital separation and period is the combined 
mass of the system. The orbital period and radius per Mo of system 
mass for the VV V-WIT-08 binary system are given in Table 4. We 
obtained an estimate of the mass of the giant in Section 3.2.1 of 
— 1 Mo for the isochrone model. If the giant has a smaller radius 
than measured from the isochrones, as in the bulge model, then this 
requires lower mass to maintain the value of log g measured from 
the NIR spectrum. The mass implied for the giant star in the bulge 
model is of order 0.4 Mo, with large error bars due in part to the 
large assumed systematic error on log g. 
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As noted in Section 4.2.3, the inferred upper limit on major 
axis radius and axis ratio are dictated largely by the prior. Only 
the minor axis radius is strongly constrained by the data, since 
the depth is almost 100 per cent and the light curve is not flat- 
bottomed the radius of the occulter on its shortest dimension must 
be approximately equal to the radius of the giant. This caveat should 
be kept in mind when considering the size of the occulter. The lower 
bound of the inferred major axis radius is constrained to areasonably 
high degree of confidence, due to the need for it to be somewhat 
larger than the minor axis radius, since occulters that are circular 
in projection cannot produce the observed light curve (as seen in 
Figure 8). Relatedly, it’s possible that we are instead seeing an edge- 
on disc and that what we had thought of as the minor axis radius is 
instead its thickness. A conservative estimate of the maximum size 
of the disc is the orbital separation of the system, which is system 
mass and giant radius dependent. 


The J band event depth as measured from the models is 
9593" per cent at lo, the maximum possible transparency 
of the occulter is therefore 3 per cent, which is the upper limit dic- 
tated by the event modelling. This is inconvenient as achromaticity 
is easier to explain with a high optical depth, but we note that swap- 
ping occulter transparency for flux contamination will yield the 
same light curve. As noted previously, flux contamination may have 
evaded our test in Section 4.3 due to the signal being of roughly 
the same magnitude as the photometric uncertainties. Addition- 
ally, the contaminating flux test did not include contaminants with 
T.g < 5000 K, leaving room for a cooler companion. Furthermore, 
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Figure 11. Posterior pdf of giant radius and occulter size and velocity 
for the two giant models considered: one obtained from isochrone fitting 
(isochrone model), and one obtained assuming a bulge location for the giant 
(bulge model). 


Table 4. Predicted occulter parameters for the two giant star scenarios we 
considered. ‘bulge model’ is a sub-luminous giant star at the distance of the 
Galactic bulge, and ‘isochrone model’ is a giant star at ~ 20 kpc with an 
improbably large space velocity. We give the median and 95 per cent con- 
fidence interval for each parameter. Subscripts ‘g’ and ‘o’ refer to the giant 
star and occulter respectively. Most of the parameters are self explanatory. 
Min Mx,,o is an approximate minimum Ks band absolute magnitude of 
the occulter under the assumptions (i) that it completely obscures the giant 
and therefore constitutes 3 per cent of the combined flux of the giant and 
occulter, and (ii) that the occulter is not itself obscured. The orbital period 
and separation axis estimates Porbit and agrpit are per Mo of system mass, 
under the assumption of a circular orbit. 


Parameter bulge model isochrone model Units 
Tg 0. 14*007 0,3220 07 au 
Tmaj,o O27 34 0.54+9:64 au 
Tmin, o 0. 14*007 0.30*005 au 
Viel,o 4.855 10:5*17 kms! 
Pori/ Mosg 242707 2394 years M5! 
&orbit/ Mog 39720 8*5 au Mg! 
min M,.o 0.34+0:83 -1 22 mag 


there is plenty of room in the event model posterior for opaque 
occulters. 


5.2 The nature of the occulter 


We have demonstrated that a solid occultation model fits the avail- 
able data. We cannot rule out an occulter with a radial transparency 
gradient, for example, but this could not explain the relative achro- 
maticity across wavelengths of 0.5 to 2.2 um unless it were com- 
prised principally of large-grained material. Further investigation, 
including modelling of the event with such an occulter would still 
be worthwhile. 

As noted in Sections 3.2.1 and 3.2.2, the measured radius 
and/or space velocity of VV V-WIT-08 are problematic to some de- 
gree. If the giant is on an ordinary single-star evolutionary pathway 
then it has an unusually large space velocity, whereas if it has a 
typical space velocity and distance for objects in this direction then 
it has an unusually small radius. It is possible there is a middle 
ground to be found between the two scenarios. An explanation for 
the event which also offers an explanation for the apparent unusual 
nature of VV V-WIT-08 is desirable. 

The required features for the occulter are that (i) it is bound 
to the giant, (ii) it is faint, (iii) it has a radius greater than 0.25 au, 
(iv) it is elliptical in projection. Additional desirable features are 
that (v) its presence also explains a high space velocity or smaller 
giant radius, (vi) it can explain a mildly blue / — Ks colour during 
the first half of the event. 

Based on these requirements, in particular requirement (iv), it 
is likely the occulter is a circumstellar disc. A circumstellar disc 
could be hosted by a body which is gravitationally bound to the 
giant, and its hosting body could be faint or itself obscured by the 
disc. Dusty circumstellar discs are observed around all phases of 
stellar evolution (Zuckerman 2001). We discuss them in turn. 


5.2.1  Pre-main-sequence disc 


We consider the possibility of occultation by a pre-main-sequence 
or YSO disc. Brown dwarf discs found in star forming regions 
typically have radii ranging from a few au to several tens of au (e.g. 
Scholz et al. 2006; Hendler et al. 2017; Rilinger et al. 2019), which 
is consistent with the inferred radius of the occulting disc in our 
model. 

One notable example is Mamajek's object, a~ 16 Myr old pre- 
main sequence star that exhibited a long, deep and complex dipping 
event which is consistent with occultation by a series of rings akin 
to the densest B-ring of Saturn (Mamajek et al. 2012). These rings 
are hosted by an unobserved companion body, likely to be a brown 
dwarf or a massive planet. 

For Mamajek's object the outermost ring is approximately 0.6 
au in radius and, assuming similar properties to Saturn's B-ring, the 
entire system must contain ~ 100 Myoon of material (Kenworthy 
& Mamajek 2015). The opacity of the rings appears to increase 
inwards and it is believed the sequence of rings was potentially 
carved out by exomoons. It is plausible that our occulter is similar 
to the disc in Mamajek's object, based on the size and the range 
of opacities, and so would be associated with a similar mass of 
material. Since for our object the light curve is smooth the occulter 
cannot be as structured as Mamajek's. 

YSO accretion discs have a typical lifetime of a only few Myr. 
We would not expect to find a bound YSO companion to the giant 
unless the giant were very massive, distant and with Vtan >> Vesc. 
Furthermore, if we adopt a bulge location for the giant star, then the 
occulter is at a heliocentric distance of at least 6 kpc and therefore 
greater than 400 pc below the plane of the Milky Way. The distance 
from the plane is naturally larger for the isochronal giant model 
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with a larger heliocentric distance. The scale height of molecular 
clouds and pre-main sequence stars in the Milky Way is only 35 to 
50 pc (e.g. Stark & Lee 2005; Bobylev & Bajkova 2016), so such a 
youthful object is very unlikely to be found in such a location. 

The difficulties with the pre-main sequence disc scenario con- 
tinue with the colour of the event. Extinction by optically thick 
YSO discs is generally not achromatic: at optical and near infrared 
wavelengths sub-micrometre-sized grains are generally observed 
(D'Alessio et al. 2006). Larger grains settle to the disc's mid-plane 
and are observed only at far-infrared and submillimetre wavelengths. 
Grey extinction by larger grains has occasionally been suggested at 
optical and near infrared wavelengths (e.g. Throop et al. 2001; 
Bozhinova et al. 2016) but to our knowledge it has never been 
confirmed by follow-up work (see e.g Miotello et al. 2012). 

On consideration of the above we conclude that our occulter 
cannot be a disc hosted by a pre-main-sequence star. 


5.2.2  Main-sequence discs 


Large and optically thick discs hosted by main sequence stars are 
generally confined to those that are relatively young, and these can 
be discounted using the same scale height argument as for pre-main- 
sequence stars. 

Older main-sequence stars can host debris discs, these were 
first observed around Vega and are analogous to the Solar Sys- 
tem's Kuiper belt. Infrared excesses suggest debris discs are present 
around many main sequence stars (Hughes et al. 2018). These discs 
extend up to 10 — 200 au from the host star, and since the upper limit 
on occulter Fmaj is poorly constrained by our data we cannot rule 
out this scenario on these grounds alone. However, main-sequence 
star debris discs are optically thin, so we consider an occultation by 
a debris disc hosted by a companion main-sequence star to be an 
unlikely scenario. 


5.2.3 White dwarf debris disc 


Farihi (2016) reviews the evidence for circumstellar discs around 
white dwarfs. The pollution of metals in white dwarf atmospheres 
suggest many of these objects experience accretion. Jura (2003) 
hypothesised the accretion could be from an optically-thick disc of 
dust formed from the tidal disruption and subsequent grinding down 
of asteroids. Such a picture requires the disc to extend no further than 
the white dwarf’s Roche radius of ~ 1 Ro, which is corroborated by 
infra-red excess measurements (e.g. Jura et al. 2007). Furthermore, 
several white dwarfs have been observed to have double-peaked Call 
emission lines indicative of a gaseous component in the surrounding 
disc constraining the outer radius to ~ 1 Ro (e.g. Gansicke et al. 
2006). 

Whilst typical white dwarf debris discs disrupted at the Roche 
limit are too small, the expansion of a red giant causes major dis- 
ruption of the associated solar system due to the loss of mass from 
the system. This might lead to increased collisions amongst the 
remaining asteroids and planets in a solar system (see modelling 
by Maldonado et al. 2020) that could perhaps cause the formation 
of a much larger white dwarf debris disc than has been observed 
hitherto, but this is mere speculation. 

From consideration of the timescale of Poynting-Robertson 
drag on small particles and the observed atmospheric pollution, it 
is understood that the discs must be optically thick (Rafikov 2011), 
although Bonsor et al. (2017) argue this may be inconsistent with 
the frequency of infrared excesses observed in an unbiased sample 
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of white dwarfs. Finally, theoretical models suggest that sublimation 
occurs around 0.2 Ro (Rafikov 2011) producing an inner hole, whilst 
the outer edge of the disc is sharply truncated as optically thin ma- 
terial is rapidly dragged into the higher density regions (Bochkarev 
& Rafikov 2011). The event modelling suggests that the impact pa- 
rameter ro is always smaller than rg, i.e. the center of the occulter 
passes over the stellar disc. However, the log-uniform prior on ro 
tends to restrain it to smaller values. It is possible there is room 
for models with large rgaj and rg. Such models could have an inner 
hole in the occulting disc which would not impact the observed light 
curve. 

Based upon these considerations, although the optical thick- 
ness, hard edge and low luminosity favour a white dwarf debris 
disc, they do not match the observed size. The conservative mini- 
mum size of the occulter is 50 Ro. We therefore conclude a debris 
disc hosted by a white dwarf is an unlikely potential occulter. 


5.2.4  Post-main-sequence — B-type subdwarf disc 


B-type subdwarf (sdB) stars are core helium burning stars at the 
blue end of the horizontal branch. Their existence can be explained 
through binary population synthesis models, around half are found 
in close binaries with white dwarfs or low-mass main-sequence stars 
(Heber 2009). 

The scenario favoured by Rodriguez et al. (2016) in the case 
of TYC 2505-672-1 — the most similar example to VV V-WIT-08 
that we have found in the literature — is of an sdB star hosting a 
large and opaque disc. In their case, UV observations by GALEX 
(Bianchi et al. 2011) facilitated the identification of the sdB star 
in the composite SED. The GALEX All-sky Imaging Survey and 
Medium-depth Imaging Survey cover V V V-WIT-08 but neither re- 
sulted in a detection at its position. Short wavelength observations 
are hampered by the relatively high extinction in this direction in 
any case. 

While we consider this a plausible scenario for VV V-WIT-08, 
our concern is that the formation of an sdB star typically requires 
mass transfer from the sdB progenitor to a companion. If this were 
the case, it would make an apparent reduction in the size of the giant 
star even harder to explain, though perhaps a suitable explanation 
can be found should further observation reveal the presence of such 
a companion. 


5.2.5 Disc comprised of matter stripped from the giant 


We speculate that mass transfer from the giant to the gravitational 
influence of a companion with sufficient angular momentum might 
explain the presence of a disc. This scenario is attractive as it offers 
a convenient explanation of the potentially smaller than expected 
radius of the giant. Any plausible companion of the giant star might 
be a host to such a disc — a main sequence star, itself a close binary 
system, a white dwarf or other compact object — but there are some 
additional considerations. 

If the companion is on the main-sequence it naturally must 
have been less massive than the giant during its main-sequence in 
order for the giant to have evolved faster. Most main-sequence stars 
less massive than the giant are fainter than our detection threshold 
even if they are not themselves obscured. We note that the canonical 
€ Aurigae scenario is of a born-again post-AGB F-type supergiant 
primary that has a B-type main sequence companion hosting a large, 
opaque disc. While the spectral types of the stars are different and the 
disc is bright at long wavelengths (Hoard et al. 2010), the e Aurigae 
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system resembles VV V-WIT-08 in many ways. An explanation for 
the lack of a mid-infrared excess in the AIIWISE and GLIMPSE 
photometry (see Figure 5) of VVV-WIT-08 might simply be that 
the disc is too small. The surface brightness of 3600 K blackbody 
(the giant star) is ~ 50 x larger than that of a 500 K blackbody in 
the Spitzer 8.0 um band, and ~ 150 x brighter in the 5.8 um band. 
With the GLIMPSE uncertainties of ~ 0.05 mag for VV V-WIT-08 
we would not expect to be able detect the presence of a 500 K 
disc at 8.0 um if it has an angular area smaller than ~ 2.5 x that 
of the giant. This not considering the potential for contamination 
in the relatively low resolution mid-infrared data due to the close 
proximity of unrelated neighbouring stars on the sky. 

Another consideration is that our estimate of the current orbital 
separation of the two components is too large to facilitate Roche 
lobe overflow. Our predicted orbital separation per Mo of system 
mass is given in Table 4. Using the posterior on the mass and 
radius of the giant from Sections 3.2.1 and 3.2.2, and the Eggleton 
(1983) approximation to the Roche lobe radius we find the closest 
the giant comes to filling its Roche lobe is ~ 10 per cent with 
Mgiant/Mocculter = 4 ~ 1 for companion masses between 107? Mo 
and 10 Mo. The more massive the companion, the smaller the filled 
fraction of the Roche lobe radius of giant star, due to the larger 
orbital separation required to explain the velocity of the occulter 
relative to the giant. Hence, unless the companion is on a highly 
elliptical orbit, outward migration due to mass loss from the system 
is probably required. 

We should also consider that the event modelling prefers an 
occulter with a position angle that is significantly offset from the 
direction of velocity (0 ~ 60°). If the disc were comprised of 
material stripped from the giant we might expect it to be aligned with 
the orbital plane of the companion. This though might be explained 
through a companion with an inclined rotation axis relative to its 
orbital plane and a magnetic field facilitating transfer of angular 
momentum to the accreted material, for example. 


5.2.6 Black hole or neutron star disc 


Related to Section 5.2.5, one possibility might be an accretion disc 
around a black hole or neutron star companion. Lasota (2016) 
estimates the size of black hole accretion discs as up to 10° 
Schwarzschild radii. A disc of minimum size 0.25 au therefore 
requires a black hole with a minimum mass of approximately 
~ 10 Mo. If we adopt the more realistic 1 au minimum occul- 
ter radius then the minimum mass increases to ~ 50 Mo. Such 
accretion discs are likely optically thick but will usually emit X- 
rays, yet no sources are present within 2’of VVV-WIT-08 in either 
the XMM-Newton Serendipitous Source Catalogue (3XMM-DR8; 
Rosen et al. 2016) or the Second ROSAT all-sky survey (2RXS) 
source catalog (Boller et al. 2016). Even if quiescent now, quiet pe- 
riods are punctuated by outbursts, so the object would nonetheless 
be a transient X-ray source. However, there are no very sensitive 
X-ray observations of this field, so it is conceivable that such an 
object might have been missed. 

We discuss in Section 5.2.5 the increase in the predicted orbital 
separation of the VV V-WIT-08 system components that is caused 
by such a massive companion, and hence the inability of the giant 
to fill its Roche lobe. This then requires a quiescent accretion disc 
that is formed of material stripped from the giant prior to a radial 
migration, or is formed of remnant material in the vicinity of the 
companion due to mass loss during its evolution to a compact object. 
The latter scenario leads us to considering discs formed from black 
hole fallback material of the type described by Perna et al. (2014) 


and elaborated on in the context of stellar occultation by Wright 
& Sigurdsson (2016). Such a disc is a plausible occulter in the 
VVV-WIT-08 system. We note that black hole fallback discs require 
an inner hole, and this is not necessarily ruled out by our occulter 
modelling. 


5.3 Contextualising VVV-WIT-08 


In addition to VVV-WIT-08 we have identified two further promis- 
ing candidate events in our VVV and VV VX data based on shallow 
searches, VVV-WIT-10 and VVV-WIT-11. These additional can- 
didates reside on the giant branch of the CMD. VVV-WIT-08 is 
a confirmed late-type giant star, the only open question being its 
distance and hence radius. The analyses of TYC 2505-672-1 by 
Rodriguez et al. (2016), and ASASSN-21co by Way et al. (2021) 
indicate that they too are late-type giant stars. While the € Auri- 
gae primary is a somewhat hotter FO star with supergiant surface 
gravity. We speculate that long-term, wide-field stellar photometric 
monitoring is beginning to reveal the presence of a population of 
long-period eclipsing binaries composed of giant stars and opaque- 
disc-hosting companions. 


6 CONCLUSIONS 


We have reported and analysed an unusual dipping event in the light 
curve of a late-type giant star, VV V-WIT-08. The dip is unusually 
long at over a few hundred days and has a depth of 97 per cent. It is 
largely achromatic, though we observe a mildly bluer than expected 
I — Ks colour during the first half of the event. 

The object was found serendipitously in the VISTA Variables 
in the Via-Lactea (V VV) survey. Subsequent searches revealed two 
further event candidates, V V V-WIT-10 and VV V-WIT-11, although 
their light curves are much less well sampled. Contreras Peña et al. 
(2017) and Lucas et al. (in preparation) find a number of long 
duration ‘dipper’ events associated with Young Stellar Objects in 
VVV searches for high amplitude variable sources, but none share 
the distinctive smooth dimming of VV V-WIT-08. 

The proper motion of VV V-WIT-08 is entirely consistent with 
a location in the Galactic bulge. However, the fitting of isochrones 
to the photometry gives a larger distance of ~ 20 kpc. Such a 
distance seems relatively unlikely, as the inferred space velocity 
is then comparable or larger than the Galactic escape speed. It is 
possible that VVV-WIT-08 is relatively underluminous, perhaps a 
sub-subgiant star or other product of binary stellar evolution, which 
renders the isochronal fit suspect. 

The behaviour of the light curve of VVV-WIT-08 does not 
correspond to any known intrinsic stellar variability. It is almost 
certainly an occultation of the giant star. We find that an occultation 
due to a chance alignment with the giant star requires an improbably 
large space density of dark foreground objects. Hence the occulter is 
most likely gravitationally bound to the giant. By modelling the light 
curve as an occultation by an object that is elliptical in projection and 
has a uniform transparency, we find that the occulter must be faint, 
likely optically thick and possess a radius or thickness in excess 
of 0.25 au. The inferred orbital parameters of the system weakly 
constrain the maximum size of the occulter to tens or hundreds of 
au, dependent on the true radius of the giant and the mass of the 
system. Our modelling implies that a hard edge of the occulter is 
necessary. We justify this on the grounds that it is more difficult 
to explain the relative achromaticity with an opacity gradient, but 
modelling of such a system is encouraged nevertheless. 
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We speculate that V V V-WIT-08, the two additional candidates, 
TYC 2505-672-1 and ASASSN-21co might be representative of a 
population of long-period eclipsing binaries composed of late-type 
giant stars and opaque-disc-hosting companions. e Aurigae might 
be considered an additional member if the spectral types of its com- 
ponents can be overlooked. If this population exists, future wide- 
field time series photometric surveys such as the Legacy Survey of 
Space and Time (LSST) at the Vera Rubin Observatory (Sajadian 
& Poleski 2019) may reveal more such objects. 

We considered a number of possible astrophysical objects as 
candidates for the occulter. Debris discs around main sequence stars 
are too optically thin. Whilst white dwarf debris discs are optically 
thick, they are too small. Accretion discs around black holes and 
neutron stars usually emit X-rays, but a black hole fallback disc of 
the type described by Perna et al. (2014) might plausibly explain 
the occultation. A disc formed of material stripped from the giant 
is an attractive scenario. A difficulty with this is that the present 
orbital separation cannot facilitate Roche-lobe overflow. However, 
this could be explained by an increase in orbital radius due to mass 
loss from the system. An optically thick YSO disc or an analogue 
of Mamajek's object are approximately the right size and could be 
optically thick. However, the height of the object above the Milky 
Way disc plane argues against both possibilities, given the scale 
height of young populations. Despite intensive efforts, it is clear 
that we have left room for further work on this intriguing object! 
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